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RADIATIVE AND GAS COOLING OF FALLING MOLTEN DROPS 

INTRODUCTION 


Molten drops of various sizes of selected materials can be cooled and 
solidified within the 2. 5 s of free-fall time available in the Marshall Space Flight 
Center ( MSFC) drop tube facility* For materials with a high melting tempera- 
ture, above approximately 1800 K, radiation cooling 13 sufficient to completely 
solidify the drops of small diameter (1 to 5 mm). However, for drops of mate- 
rials with lower melting points, heat loss by radiation cooling alone will not 
allow complete solidification within the 2. 5 s time limit. Therefore, it becomes 
necessary to back-fill the drop tube with known pressures of a cooling gas to 
increase the cooling rate of a falling drop. Helium gas would seem to be a good 
candidate since this gas is inert, readily obtained, and ims a high vr’ue of 
thermal conductivity relative to other gases. The purpose of this report is to 
provide a good estimate of the effect that the presence of helium gas in the drop 
tube would have upon the rate at which heat will be lost by a falling drop. This 
effect should be evidenced by a change in the time required to solidify a falling 
drop and by a change in the supercooling rate of the molten drop if the drop does 
not begin to solidify after release. 


THEORY 


The rate at which a heated sphere will lose heat energy during free fall 
in the drop tube filled with a gas is given by Zemanskv [1] as: 


^ = -eA<r(T< -T 4 ) -hA(T -T ) . (l) 

dt o o 

The first term on the right of equation (1) is due to radiation loss and Is depen- 
dent on the emissivity of the material, e; the surface area of the sphere. A; the 
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Soltzman' s constant, a ; as well as the fourth power of the sphere temperature, 
T, and the ambient temperature, T q . The second term on the right of equa- 
tion (l) is due to a combination of conduction and forced convection and is 
dependent upon A, T, T^, and the heat transfer coefficient, h. The amount of 

heat a drop must lose for complete solidification will be equal to the product of 
the mass, m y of the drop and the latent heat of fusion, H^, of the drop material. 

If during solidification the temperature of the drop is assumed to be constant 
at the melting point, T , of the drop material, then equation (l) can be reduced 
to 


t = miy leA<r(T m 4 - T q 4 ) + hA(T - Tj\ 


( 2 ) 


The solidification time, t , of equation (2) will be the time required to completely 

s 

solidify the drop. If the drop tube is evacuated, h goes to zero and the solidifica- 
tion time becomes 




(3) 


A comparison of equations (2) and (3) for a given drop shows the effect Jiat the 
presence of a gas environment will have upon solidification time for a falling 
drop. 


If the drop does not begin to solidify after release, it will supercool. In 
this event, Q of equation (l) will be equal to the product of the mass, m, of the 
drop and the specific heat, C, of the molten material. Equation (l) will then 
reduce to the differential equation 


HT 

“T” = -K (T 4 - T 4 ) - H(T - T ) 
dt o o o 


(4) 
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where 



A notation to equation (4) has boon provided by Kat* and Wills (2) by assuming 
that the drop material parameters do not change during the time of free-fall, 
that the drop Is spherical, and that the drop Is small enough so that there Is no 
thermal lag between the sphere’ s surface and center. The solution to equa- 
tion (4) is then 






X 


t M it 


where 


K » 20 + 27H /T 
1 o o 
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and 


(5c) 


K 2 = '7k, 2 + 32 1 . 


The time t of equation (5a) is the time required for a drop to supercool from 
c 

the initial temperature, T^» ■> the final temperature, The constant is 

evaluated at t = 0 and T = T^. If the drop tube is evacuated, h goes to zero and 
equation (5) reduces to 


t 

c 


1 

4K T 3 


o o 



( 6 ) 


where the constant is again evaluated at t = 0 and T = T^. The effect that the 

presence of helium will have upon the supercooling rates of falling drops can be 
seen from a comparison of the supercooling times of equatio:. (5) and (6). 


ASSUMPTIONS 


To calculate the cooling and solidification times for a drop, certain 
assumptions must be made* As previously stated, the drop is assumed to be 
spherical, have no thermal lag between surface and center, and have constant 
material parameters during the time of cooling. The validity of assuming the 
drop to be isothermal during cooling is addressed in the Appendix. 

An empirically determined heat transfer coefficient, h, for a sphere in 
a gas flow with a Reynolds number between 17 and 17 000 is recommended by 
McAdams [3 ] as 
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(7) 


where D is the diameter of the drop, V Is the free stream velocity of the gas 

CD 

flow, and and are the thermal conductivity and kinematic viscosity of the 

gas at the thin ftlm temperature, T . Although the drop is accelerated through 

K 

the gas by gravity when tolling in the drop tube, for ease of calculation, the 
velocity of the drop relative to the gas will be assumed to be a constant 12 m/s, 
which is the average velocity of a drop falling the length of the drop tube 
(30.5 m). Although, according to Loeb (4), the thermal conductivity and 
viscosity of the gas are pressure independent for the pressures to be considered 
in this report (i. e. , 1 to 760 Torr), those parameters are evaluated at T_ 

n 

since they are temperature dependent. Since equation (7) is dependent on the 

kinetic viscosity which is the quotient of the viscosity and density of the gas, the 

heat transfer coefficient h will be pressure depends f . Holman [5j gives the 

thin film temperature as the average of the surface ; mperature, T , of toe 

8 

sphere and the free stream temperature. T . of the gas 

T + T 
a • 
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Since the change with temperature of K and V is slow, T can be assumed to 

t ( 8 

be the average of T and and equation ($) becomes 
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T, + T . + 2T 
1 f o 

9 


(9) 


bv also assuming T = T . 

® o 

Although the previously described assumptions Introduce some error Into 
equations (2), (3), (5), and (t>), these equations should s.tll provide a good 
estimate of the supercooling and solidification times for a drop In both a helium 
and vacuum environment. 
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RESULTS 


Niobium, copper, and lead have been selected as typV'al materials 
because of the wide range of material parameters, especially the range of 
melting points, characteristic of these materials. The range of melting points 
of these materials is very important cor jiderlng the fact that the drop is at or 
slightly above the melting point when the drop is released. The characteristic 
parameters for the selected materials are presented in Table 1. 

The solidification times, as given by equation (3) for radi°tlon cooling 

only, for niobium, copper, and lead drops of 3, 5, 7, and 10 mm diameter are 

present*! in Table 2. The solidification times of Table 2 reflect the fourth 

power ot T dependency of t of equation ( 3) as evidenced bv the marked 
ms 

increase of the required times for complete solidification for the drops of ne 
materials with the lower melting points. Also from Table 2 it should be noted 
that radiation cooling is not sufficient to solidify drops of copper or lead within 
the 2. 5 s time limit of the drop tube. 

The solidification times for niobium, copper, and lead drops falling 

through a helium gas, as given by equation (2), are presented in Figures 1, 2, 

and 3 as a function of helium gas pressure. The helium pressure dependency 

of the solidification time enters equation (2) by way of the kinematic viscosity 

of the gas from the heat transfer coefficient. As can bo seen in Figure 1, the 

effect of helium cooling on the solidification times for niobium drops would be 

sufficient to allow solidification of 7 and 10 mm diameter niobium drops in the 

dt\ tube. The cooling effect of helium is even more pronounced in the curves 

of Figures 2 and 3, which show that copper drops of up to 5 mm in diameter and 

lead drops of diameters up to 7 mm could bo completely solidified if the drop 

tube were filled with helium gas at near standard atmospheric pressure. It 

should be noted from Table 2 that none of the copper and lead drops discussed 

would solidify in the drop tube if it were evacuated, A comparison of Figures 1, 

2, and 3 reveals that the cooling of the lower melting point copper end lead drops 

is less dependent on radiation than on heat transfer to the helium, as evidenced 

by the rapid decrease in drop solidification times as helium gas pressure 

increases. To illustrate the fact further, the percenta.ee of heat loss at T due to 

m 

the radiation term of equation ( 1) is presented in 1 igure 4 fi..- the three materials 
as a function of helium pressure. Fig . e 4 clearly shows that radiation cooling 
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Is much less a factor for lead and copper drops than for the higher melting point 
of niobium drops and, therefore, agrees with Figures 1, 2, and 3, Since radia- 
tion loss adds very little to the cooling of copper and lead drops near the stand- 
ard atmospheric pressure of helium, the solidification times for these drops 
should depend less on the melting temperature than on other material parameters. 
Therefore, the fact that larger lead drops can be solidified in the drop tube than 
copper drops (as can be seen by comparing Figures 2 and 3) can be attributed 
to the relatively low heat of fusion of lead. 

It should be noted that the heat transfer coefficient used In the calcula- 
tions presented ir. Figures 1 through 4 is not accurate for Reynolds numbers 
below 17. The regions of the curves in Figures 1 through 4 for which the 
pressure of the helium gas is such that the Reynolds number of the gas flow past 
the foiling drop is near or below 17 are shown as dotted. These regions are 
presented since they should still provide a good estimate of the effect the helium 
will have upon drop solidification times for this range of pressures. 



Also of interest is the rate at which a drop will supercool if it doe* not 
begin to solidify immediately after being released. The times required for 
3, 5, 7, and 10 mm diameter drops of niobium to supercool a given amount are 
given in Figures 5 through 8. These are times for niobium drops falling through 
various pressures of helium gas, as given by equation (5), and times for drops 
falling through a vacuum, as given by equation (6). Corresponding supercooling 
curves for copper are given in Figures 9 through 12, and for lead in Figures 13 
through 16. Although the drops are not expected to supercool more than a few 
hundred degrees, the curves of Figures 5 thi ough 12 have been extended to 
1000 deg of supercooling to further illustrate the effect of the presence of helium 
gas on the supercooling rates of niobium and copper drops. 



From Figures 6 through 9 it can be seen that the niobium drops will 
supercool very quickly. Before a .’> mm niobium drop has fallen 1 m (0.45 s) 
in vacuum, it could have supercooled as much as 250 deg. Within the 2.5 s limit 
of the drop tube, a 3 mm diameter drop of niobium can possibly supercool down 
to the homogeneous nucleation tomperature which is given by Miller ( 8] as 
approximately 700 K. Helium gas cooling would increase thv supercooling rate 
by as much as 80 percent near the atmospheric pressure of helium. This 
Increased cooling rate could also be expected for the 5, 7, and 10 mm diameter 
niobium drops. For the copper drops discussed in Figures 10 through 13, the 
supercooling rates are much slower. A 3 mm diameter copper drop 
possibly supercool approximately 50 dog in the evacuated drop tube. However, 
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with helium cooling, the same size drop could supercool to the homogeneous J 

nucleatlon temperature of 240 K alter only 0. 54 s of free fail. Figures IS j j 

through 16 show an even greater difference between the vacuum and helium j 

curves for lead drops. If the drop t*»be were evacuated, a 3 mm diameter lead I 

drop could not be expected to supercool more than a few degrees. Likewise, • 

with helium present In the drop tube at standard atmospheric pressure, the 
same diameter lead drop could possibly supercool to the homogeneous nucleatlon 
temperature of 80 K after only 0. 34 s of free fall. By comparing the super- 
cooling curves of niobium, copper, and lead, it can be seen that cooling Is more t 

dependent upon heat transfer to the helium gas than on radiation loss for the j 4 

lower melting temperature materials. This result could have beon predicted 
from Figure 4. 

CONCLUSION 


Since drops of samnle materials are released In the drop tube at or 
slightly above their melting temperature, the effect of helium gas cooling will 
be greater for drops of materials with low melting points. By the use of helium 
of less than 1 atm of pressure, solidified drops of significant size (l.e. , 

> 3 mm) can be obtained for materials which could not be solidified in useful 
sizes in the 2. 5 s time limit for the evacuated tube. Also, if they do not begin 
to solidify immediately upon release, the rate at which the molten drops super- 
cool will be considerably afTected. The presence of helium can increase the 
supercooling rale enough to supercool drops to the homogeneous nucleatlon point. 
Thus, the present drop tube facility can be used to study nucleatlon and the 
solidified strut tore of drops of metals and alloys supercooled the maximum 
extent possible. In conclusion, it can be stated from the foregoing estimates 
that the presence of helium gas in the drop tube would have a significant effect 
on both time required to solidify drops of sample materials and on the super- 
cooling rates of the drops. Thus, the use of helium gas could greatly extend 
the versatility and usefulness of the MSFC space processing drop tube to carry 
out containerless low-gravity supercooling and solidification experiments. 
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TABLE 1. CHARACTERISTIC PHYSICAL PROPERTIES 
OF SAMPLE MATERIALS (6,7) 



Niobium 

Copper 

Lead 

Density 
(g/ cm 3 ) 

8.60 

j.96 

11.34 

Melting Point 
(K) 

2741 

1357 

600 

Heat Capacity 
(J/g-K) 

0.268 

0.385 

0.126 

Heat of Fusion 
(J/g) 

284.6 

211.8 

26.4 

Emissivtty 

(NearT ) 

m 

0.25 

0.16 

0.075 


TABLE 2. SOLIDIFICATION TIMES IN SECONDS FOR DROPS 
FAMING IN A VACUUM 


Dtametor 
( mm) 


Niobium 


Coppor 


ao.9 


loa.o 


Lead 


288. 3 
480.-1 
072.6 
960.9 
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Figure 1. 
of 3, 


Time required to completely solidify niobium drops 
, 7, and 10 mm diameter as a function of helium 
gas pressure. 
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Figure 3. Time required to completely solidify lead drops 
of a, 5, 7, and 10 mm diameter as a function of helium 
gas pressure. 
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Figure 5. The amount of supercooling for a 3 mm diameter niobium 
sphere plotted as a function of time before nucleatlon. (The 
calculations are for a molten drop falling in a vacuum 
or helium gas at different pressures. ) 
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SUPERCOOLING (K) 


Figure 6, The amount of supercooling for a 5 mm diameter niobium 
sphere plotted us a (Unction of time before nucloatlon. (The 
calculations are for a molten drop falling In a vacuum 
or helium gas at different pressures.) 




SUPERCOOLING (K) 


Figure 7. The amount of supercooling for a 7 mm diameter niobium 
sphere plotted as a function of time before nucleatlon. ( The 
calculations are for a molten drop foiling In a vacuum 
or helium gas at different pressures. ) 





SUPERCOOLING <K) 


Figure 8. The amount of supercooling for a 10 mm diameter niobium 
sphere plotted as a function of time before nucleatlon. (The 
calculations are for a molten drop falling In a vacuum 
or helium gas at different pressures.) 
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Figure 9. The amount of supercooling for a S mm diameter copper 
sphere plotted as a (Unction of time before nucloation. (The 
calculations are for a molten drop falling in a vacuum 
or helium gas at different pressures.) 
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Figure 10. The amount of supercooling for a 5 mm diameter copper 
sphere plotted as a function of time before nucleatlon. (The 
calculations are for a molten drop falling In a vacuum 
or helium gas at different pressures. ) 
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7 mm COPPER 
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SUPERCOOLING IK) 

Figure 12. The amount of supercooling for a 10 mm diameter copper 
sphere plotted as a function of time before nucleatlon. (The 
calculations are for a molten drop falling in a vacuum 
or helium gas at different pressures.) 
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3 mm LEAD 



SUPERCOOLING (K) 

Figure 13. The amount of supercooling for a 3 mm diameter lead 
sphere plotted as a function of time before nucleation. (The 
calculations are for a moltei: drop falling in a vacuum 
or helium gas at different pressures.) 
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ralrulationa aro for a motion drop falling In a vacuum 
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SUPERCOOLING (KJ 


Figure IS. The amount of supercooling for a 7 mm diameter lead 
sphere plotted as a (Unction of time before nucleatlon. ( The 
calculations are for a molten drop falling in a vacuum 
or helium gas at different pressures.) 
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Figure 16, The amount of supercooling for a 10 mm diameter lead 
sphere plotted as a function of lime before nucleatlon. (The 
calculations are foi a molten drop falling in a vacuum 
or helium gas at different pressures.) 
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APPENDIX 


VALIDITY OF ISOTHERMAL COOLING 


The question of the validity of the isothermal approximation, for the 
purpose of this report, can be answered by considering the temperature distri- 
bution within the cooling sphere. The radial heat flow within the sphere must 
satisfy Poisson's equation [ 9 ]. 


8*T t 2 aT a dr 
2r* r 3r ~ dr 


(A-l) 


where tr = K/ pC is the thermal diffusivity of the material and K denotes the 
material's thermal conductivity and p the material's density. The boundary 
conditions would be: 


T = Tj at t = 0 


and at the surface 


?T 
dr " 


h(T-T )+ €o(T*-T<) 
o o_ 

K 


for t>0. 


(A-2) 


(A-3) 


Since equation (A-3) is a time -dependent, nonlinear boundary condition, 
equation (A-l) has no closed form solution. However, equation (A-3) can be 
linearized by using a Taylor's series expansion which yields 


? T _ 
? r 


h(T.-T ) + ea (T, 4 -T 4 ) 
1 o » o 

K 


h + 4€oT* 
K 


(T-T,) 


(A-4) 


where T, is a temperature selected near the mid-range between T. and T^. Tj 

can be chosen so as to minimize the difference between equations (A-3) and 
(A-l) over the desired temperature range. Using boundary conditions (A-2) 
and (A -4), Carslaw and Jaeger [ 9l offer a series solution to equation (A-l) 
as: 
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^[rf^ + aJtaJ-i)]] 8111 <a V 8ta <r V 


(A-5) 


where a is the radius of the sphere and 


v = T - Tj + Tj 


V = T - Tj + Tj 


with 


T 2 = 


h(T. -T )+ Cct( T 4 -T 4 ) 
1 o * o 


KJ 


and 

, h + 4€oT, 5 

j= 

and where 0 N is given as the roots of 
a$ N cot (a$ N ) + aJ - 1 =0 . 


The temperature difference between the drop's surface and center, 
defined as AT, can now be found by using equation (A-5) and the values of 
K [lO] and a as given in Table A-l. The results show that 3 mm diameter 
copper and lead drops w ill have a AT of less than 4 K even with the drop tube 
filled to one atmospheric pressure of helium. These AT's for copper and lead 
correspond to less than 1 percent of the melting temperature. For a 3 mm 
niobium drop falling through one atmospheric pressure of helium, the expected 
AT would be less than 45 K, or approximately 1. 6 percent of the melting tem > 
perature. A niobium drop of the same size cooled by radiation would have a 
AT of only 14 K. For convenience, the ratios of the maximum AT to the 
melting temperature, T , for 3 mm diameter drops of the three materials, 

cooled with nelium gas, are also listed in Table A-l. 

Since the temperature difference from center to surface for drops of 
niobium, copper, and lead is such a small percent of the melting temperature, 
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K(J/S-ca-K) 

0.52 

3.42 

0.34 

a(cm*/S) 

0. 11 

0. 99 

0.24 

AT max 

0.016 

0.002 

0. 007 

TM 
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